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Abstract. A systematic derivation of Boltzmann equation is presented in the 
framework of closed-time-path formalism. Introducing a new type of probe, the 
expectation value of number operator is calculated as a functional of source. Then 
solving for the source by inverting the relation, the equation of motion for number is 
obtained when the source is removed, and it turns out to be the Boltzmann equation. 
The inversion formula is used in the course of derivation. 



In this letter, we present a new approach to derive the Boltzmann equation (BE). 
There have been some works on this subject both in the framework of closed-time- 
path (CTP) formalism, || and in the framework of thermo-field dynamics. || 
These approaches have an advantage that the time dependence of the number is not 
introduced by hand. Instead, a counter-term is first introduced into the CTP or 
thermo-field Lagrangian and bare propagator is calculated. Then to determine the 
counter-term, some condition, such as the cancelation of on-shell part of the self- 
energy |l], |j| or the cancelation of pinch-singularity, Q is adopted, which leads to BE. 
Since their primal purpose is to construct the non-equilibrium perturbation theory, 
BE appears as a byproduct. It is preferable if we can derive BE more directly as an 
equation of motion (EoM) of expectation value of the number operator. Moreover in 
their approaches, the condition to determine the counter-term is of course not unique, 
and the approximation made is not so clear. In the following, a more direct approach 
is studied based on the inversion method. j| |j| A new type of probe dictated from 
the counter-term approach is introduced in course of derivation. 

Let us briefly describe the inversion method which is a systematic procedure 
to derive EoM in CTP formalism. In CTP formalism, |], Q we introduce a time 
dependent source J to probe some operator of interest, say Q{(p), which is a function 
of the dynamical variable (p. Then with the Hamiltonian H of ip, the CTP generating 
functional is defined as 

e iwwi = Tr Te"^'i Pdt( ^" Jl(t)<5) p f e*C dt ("~ J ^) (1) 
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where p is the initial distribution and T and T are the time ordering and anti-ordering 
operators, respectively. The last equality is due to path-integral representation, where 
ipi and ip>2 are respectively introduced as integral variables along the forward and 
backward time branches. 



For convenience of later discussion, let us introduce "physical" representation Q 
through Jq = |(Ji + J2), and Ja = Ji — J%- Then J a = is physical and Jq plays the 
role of external force. The expectation value of Q at time t under physical external 
source J c = J can be calculated as 



Q(t) 



SW[J A ,Jc] 
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Jc=J 
Ja=0 



Q(t) 
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This gives us the expectation value Q as a functional of external source J. 

In order to obtain the EoM of Q, we solve the relation (|^) inversely to express J 
as a functional of Q. Then setting the external source J = 0, the obtained relation 
gives EoM of Q. (Inversion method Q) Formally, the general expression of EoM can 
be written with the Legendre transformation of W. But practically, the process of 
Legendre transformation is unnecessary, and in this letter, this inversion is carried out 
in the following perturbative fashion. 

Usually Q as a functional of J is obtained as some perturbation series 



Q{t) = f[t-J] = Y^e n f {n) [t-J] 



(4) 



where e is a small parameter and f[t;J] expresses that / is a function of t and 
functional of J. Then if we write the inverted relation as 



J(t)=g[t;Q}=J2^ m 9 (m) [f,Q], 



(5) 



following simple identity is obtained 
Q(t) = f[t;g[Q}] 



= / (0) 
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where, e.g. 5f^[t; J]/SJ(s) evaluated at J = g^[Q] is abbreviated as 
Sf^Hty/fig^H 3 )- Comparing the lhs and rhs in each order of e, we obtain the 
expressions for g^ in terms of f^ n \ which we call the "inversion formulae". § 
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First of all, to make this method work, wc need a non-trivial lowest-order functional 
expression f^[t; J] which can be inversely solved for J. This becomes the key-point 
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for deriving BE. If we naively apply this method to number operator, the expectation 
value does not reveal such non-trivial dependence on J. 

Let us see the problem more closely. We consider a non-relativistic Boson 
field of a homogeneous system described by the Hamiltonian H = Hq + _ff; nt with 
#o = Ek e k^k^k' and Hi *t = lEk,k',q^k +q 4'-q^k4" where A is a coupling 
constant, which is assumed to be small and plays the role of e in (Q). Extension to 
other type of interaction is straightforward. For the initial density matrix p, we assume 
that no initial correlation exists among the different wave-number components and p 
can be written as a product form pj^., where is a density for each wave- number 

which gives n-^ti) = Trpj^l I 0j c . 

In order to derive EoM of the expectation value of the number n^{t) — $t.{t)'ip\ i (t) , 

a naive choice of the source is H — -^kW^kW^kW" Then in path-integral 
representation of CTP generating functional, this source can be built into the free 
part of the Lagrangian as 

l 7 OM - l o 7 0a 2 ) = ]T 'Wk'Yk (io) 

k 

with the matrix 

p. ( t d f ) = ( m ~ £k + JkW ° ^ (11) 

V ^ at >-\ o -m t + e k - j k (t) J ■ [LL) 



The bare propagator is essentially the inverse of the matrix in (|ll|), and with this 
propagator, if we evaluate the expectation value n^(t) in the absence of interaction, 
we just obtain the initial value {n^(t))j — nu(t{), due to the conservation of n-^ for 
A = even when J^. ^ 0. Since no dependence on J appears, we fail to obtain the 
inversion in lowest-order, and hence the inversion formulae can not be used in this 
case. Probe of the form (|TT| ) is not enough to handle the number operator. 

Then why does the counter-term method work? According to reference jjj, the 
time-local counter-term is constructed so as to keep the following structure of the full 
propagator in CTP formalism (We suppress the index of wavenumber for a while.) 



G(t, s) = -Tip 



where 'c' means the connected part and 

h(t, s) = -$(t)ft(s)) c , Ht, s) = -(^(s)m) c . (13) 

Then it turns out that the counter-term ?p*M.ijipj with the matrix 

_ / hAuQk) - ia(t) -i(ft 7 (i) - a(t)) \ 
J l{ > ~ \ i(ft 7 (t) + a{t)) -hAuo{t) - ia(t) J ' 1 j 

is allowed to be subtracted from the free part of the Lagrangian, where Au, a 
and 7 are all real functions which are determined by appropriate conditions. The 
bare propagator calculated from Lo(ipi) — £0(^2) — iPt-M-ijtpj leads to non-trivial 
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time dependence of the number in the absence of interaction. The existence of the 
parameters in non-diagonal elements is a crucial point. 

Comparing (Hm with (O), the parameter we have utilized as a physical external 



source in ( |ll[ ) corresponds to Alo in (|14J). Equation (14), however, suggests that 
another physical source corresponding to a or 7 can be introduced as a probe. Our 
choice here is the source corresponding to a. (The source corresponding to 7 can be 
treated similarly.) Then the free part of the Lagrangian including the source now has 
the matrix 

- U (t) -m t + e + u(t) )■ (15) 

Note that although the source J is introduced in this way, what we calculate in the 
following is just the expectation value of the number; We integrate ^>*(t + e)ijji(t) 
under the existence of the probe (|l5|) . (Of course other choices, e.g. ^{t — e)ip2(t), 
produce the same results.) 

From the matrix (|l5|) , the bare propagator Go is calculated by 

V(t, d t )G (t, s) = G (t, s)V(s, -d s ) (16) 
= -ihS(t-s). (17) 

Since V has been chosen as to keep the structure (12) unchanged, Go has the same 
structure in which h and k are replaced by ho and fco, respectively. Then (17) leads 
to the equations 

(ihd t -e)h (t,s)=0, (18) 

(in&-e)*o(i,«) =0, (19) 

for t > s, and 

{ihd t -e + U(t))k* (s,t) =U(t)h* {s,t), (20) 

(ihd t -e-U(t))h*{s,t)= -U(t)k* (s,t), (21) 

for s > t. The boundary conditions at t = s are given as 

h (s,s) - ko(s, s) = -1, k (s,s) - ho(s,s) = 1, (22) 
h (s,s) - hl(s,s) = 0, k (s,s) - ko(s,s) = 0. (23) 

From (23) ho(s,s) and fco(s,s) are real functions. Then the two conditions in (22) 
are identical and simply express the fact that the expectation value of the equal-time 
commutator is unity. Note that from the definition (|l3"|), ko(t, t) gives the 

expectation value of the number operator (multiplied by —1) in the absence of the 
interaction, which we denote as n^°\t). 
From (18) and (19), we obtain for t > s 

k (t,s) =e-^ ( *- s) fc (s,s) ^-/'We-i*" 8 ' (24) 
h (t,s) =e-* e(t - s) h (s,s) = ~(n^(s) + l)e-* e{t ~ s) . (25) 

Then exchanging t and s in (24) and (25) and taking the complex conjugation, h^s, t) 
and fcg(s,t) are obtained for s > t. Substituting them into (20) or (21), both equations 
turn out to give an identical result, and we find that n 1 -- ^ must satisfy the condition 

J(t) = hd t n^(t), (26) 
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which gives EoM for and is integrated as 



t;J] 



(ti) 



J(s) 



(27) 



Equations (24), (25) and (j27|) together with the structure like (12) determine the bare 
propagator. 

As already seen from (^7|) or (^6j), we succeeded to make the expectation value of 
the number depend on J in the lowest order, i.e. 0(A°). This makes the inversion 
formula applicable. The rhs of equation ( p7| ) corresponds to the desired lowest-order 
functional in (||), and (26) is the inverted relation, the rhs of which corresponds 
to g 1 - " 1 of (||). So our next task is to calculate the perturbative correction to n, and 
then to derive the correction to the EoM (f26j) with the aid of the inversion formulae. 




Figure 1. Diagram for 0(A 2 )-correction to n^.. 

With the propagator calculated above, the non-zero perturbative correction to 
nj c (i) first comes from a diagram shown in figure 0, which is of 0(A 2 ). The 
contributions of 0(A) from a tadpole type self-energy insertion vanishes due to the 
cancelation of terms from the vertices on forward and backward time branches, and 
similarly for the two tadpoles inserted diagram of 0(A 2 ). As the result 



k [t, J] = n£\t) + £ A' A*' cos (u; kiqiq ,(t' - a') 



qq 

{ ("k 0) + ^qVk + V« } - 4%VkK } + + 1) } (A (28) 



k Jy q+q'-k 
where „ „, = |(eq 



w k qq' = ll e q + e q' — e q+q'-k — e k)- R cca U that all ^ k s are functionals of 
J k given in (§7|). Equation (28) corresponds to f^+ef^+e 2 /^ of (g), where /W 
vanishes as mentioned above. 

Applying the inversion formulae, we obtain the correction to (pa) as 



A f ( ^ 

J k (i) = hd t n k (t) - — 2J / dscos (w kqq ,(f - s)^ 
q q , Jti 

x | (n k + l)(n q+q ,_ k + l)nqnq- - n k n q+q ,_ k (nq + !)( n q' + : ) } ( s )-( 29 ) 

Note that, in course of the inversion, all the functionals of J are evaluated at J k = ?iri k 
and nj^ [t; J] contained therein becomes n-^(t). If we set the external source J = 0, 
EoM for the number is obtained. The correction term is similar to the collision 
terms of Boltzmann equation, but it has non-Markovian form and contain energy 
non-conserving process. 



() 



The ordinary Markovian BE is obtained by the adiabatic expansion. Setting the 
initial time t\ = — oc, we abbreviate the products of n and n + 1 in the integrand of 
(29) as 7V (2) and expand it around the time t as N (2 ~>(s) = N ( - 2 \t) + (s-t)N ( - 2 \t) + - ■ 
regarding the time differentiations to be small. Then the integral becomes 



dscosu(t - s)N {2 \s) = irS{cj)N {2 \t) + ^N {2) (t) + ■■■. (30) 



The second term is proportional to n and gives a perturbative correction to the 
coefficient of the first term in rhs of (29), which can be neglected. Regarding all 
higher time derivatives to be small, we take into account up to the first term of 
and obtain the ordinary time-local BE with energy conserving process 



hd t n k (t) = tt\ 2 J2 S ( e q + £ q' - e q+q<-k _ e k 



q,q' 

x I (n k + l)(n q+q ,_ k + l)n q n q - - ™k™q+q'-k( n q + ^("q' + 

The key-point of our derivation is the new type of probe introduced in ( |l5| ) . After 
that, the application of inversion formula is straightforward. Of course with the usage 
of higher-order inversion formulae, we can calculate higher-order corrections to 
EoM quite systematically. This will be presented in other place. 

Physical contents of (15h becomes somewhat clear if we consider effective action 
of ip. From the CTP generating functional W with ip itself as the order parameter 
Q, the effective action T [tpA, ] is calculated through the Legendre transformation 
of W, where ipA = SW/SJc and ipe = 5W/5J&. Roughly speaking, -0A = ipi — 1P2, 
%p c — + ipz) and V is the tree part of second derivative of V. Then the source of 



the form ( |15[ ) couples to ip^tpA and corresponds to the quantity 5 2 T/Sip A SipA which is 
the 1-particle-irreducible amputated part of the correlation function {{ffi This 
may be the reason why we can handle the number with this source. Another choice of 
the source corresponding to 7 in ( |l4| ) also produces non-trivial time dependence in the 
lowest order and EoM can be derived. Although the result has somewhat complicated 
expression, it agrees with (31) after the adiabatic expansion. 

The author is grateful to professor R. Fukuda for helpful discussions. 
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